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Abstract 

We show that the classical Cauchy problem for the incompressible 3d Navier-Stokes equations 
with (— l)-homogeneous initial data has a global scale-invariant solution which is smooth for 
positive times. Our main technical tools are local-in-space regularity estimates near the 
initial time, which are of independent interest. 

1. Introduction 

We consider the classical Cauchy problem for the incompressible Navier-Stokes equation 
in R 3 x (0, oo) 

u t + uVu + Vp - Au = ) . D3 f s. 

divu = j mi? 3 x(0,oo), (1.1) 

u\ t=0 = u ini? 3 . (1.2) 

We recall that the problem is invariant under the scaling 

u(x,t) — > u\(x,t) = \u(Xx,X 2 t), 

p{x,t) p x {x,t) = X 2 p(\x,X 2 t), (1.3) 
uq(x) uq X {x) = Xu (\x) , 

where A > 0. We say that a solution u is scale-invariant if u\ = u and p\= p for each A > 0. 
Similarly, we say that an initial condition uq is scale- invariant, if uqx = uq for each A > 0. 
This is of course the same as requiring that uq be (—1)— homogeneous. 
One of our goals in this paper is to give a proof of the following result. 

Theorem 1.1. Assume u is scale-invariant and locally Holder continuous in R 3 \ {0} with 
div u = in R 3 . Then the Cauchy problem U.l\) , U.ty) has at least one scale-invariant 
solution u which is smooth in R 3 x (0, oo) and locally Holder continuous in R 3 x [0, oo) \ 
{(0,0)}. 

Previously this result has been known only under suitable smallness conditions on uq, see 
for example [5l[T5]. For small uq one can also prove uniqueness (in suitable function classes). 
It is quite conceivable that uniqueness may fail for large data. We will comment on this 
point in more detail below. 

The second important theme of our paper can be perhaps called local-in-space regularity 
estimates near the initial time t = 0. It is known that if uq G L q for q > 3, then the 
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initial value problem (II. ip . ( II. 2p has a unique local- in-time "mild solution" defined on some 
time interval (0, T), which is smooth in R 3 x (0,T) and has in many respects the same 
regularity as the solution of the heat equation in R 3 x [0,T) for times close to t — 0, see for 
example [6j[T0]. A natural question is under which condition this result can be localized in 
space: if «o is a quite general initial condition for which a generalized suitable weak solution 
u in the sense of [19] is defined and u q \b t G L q (B r ) for some q > 3, say, can we conclude 
that u is regular in £>r x [0,ti) for some time ti > 0? We prove that this is indeed the case 
under quite general assumptions, which include u which is in Lf oc and J B \u \ 2 dx —> 
for x oo. Due to non-local effects of the pressure the solution u in Bt. x [0,ti) may not 
have the same amount of regularity as the solution of the heat equation in this situation, but 
the non-local effects are limited to the influence of the "harmonic part of the pressure" in a 
suitable pressure decomposition. We can formulate this type of results somewhat loosely in 
the following statement. 

(S) Modulo the usual (and quite mild) non-local influences of the pressure, local regularity 
of the initial data propagates for at least a short time. 

We refer the reader to Section 3 for precise statements. Statement (S), in addition to 
being of independent interest, is one of the main ingredients of our proof of Theorem 11.11 

Results in the direction of (S) can be found already in the classical paper [3J. More 
recently, related questions about vorticity propagation have been studied in [27]. Our main 
result concerning (S), Theorem I3.1[ takes a somewhat different angle on (S). 

Our proof of (S) (see also Theorem 13 ,ip is based on a combination of techniques from [T2J 
[16j[T8l[T9l|25] • Heuristically, the main point is that one can obtain a sufficient control of the 
energy flux into "good regions" from the rest of the space, see Section 3. Once we know 
that only small amount of energy can move into the "good region" one can use (a slight 
modification of) partial regularity schemes in [T5JIIB] to prove regularity. 

To prove Theorem 1.1, we seek the solution u(x,t) in the form 

«(,,i) = ^ I /(i). (1.4) 
The Navier-Stokes equation for u gives 

-AU-hj--xVU + UVU + VP = 0, divt/ = 0, (1.5) 

in R 3 . For a scale- invariant u the problem of finding a scale- invariant solution of the 
Cauchy problem (II. ip . (11.21) is equivalent to the problem of finding a solution of (II. 5p with 
the asymptotics 

\U(x) — u (x)\ = o ^— ^ , x — > oo . (1.6) 

The problem ( 11. 5ft . ( 11. 6ft is reminiscent of the classical Leray's problem of finding steady- 
state solution of the Navier-Stokes equation in a bounded domain, with given boundary 
conditions for U. We will show that one can solve this problem using the Leray-Schauder 
degree theory, just as in the case of the bounded domain. The non-trivial part is to find 
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the right functional-analytic setup and establish the necessary a-priori estimates. The main 
difficulty is to find good estimates near oo. This difficulty will be overcome by applying 
statement (S) above. Heuristically it is clear that when u is given by ( II. 4p . then estimates 
of u near t = are closely related to estimates of U near oo. In Section 4 we will make this 
more precise. 

As in the case of the bounded domains, the Leray-Schauder approach gives existence of 
the solutions, but not uniqueness. In the case of bounded domains one does not generically 
expect uniqueness for large data, and this non-uniqueness is in fact expected to be quite 
typical in the context of the steady Navier-Stokes, once the data is large. Could this also 
be the case for the problem (11.51) . (11.61) ? This would lead to non- uniqueness for the Cauchy 
problem (11.11) . (11.21) with scale-invariant Mo, and by a suitable truncation of Uq at large 
\x\ possibly also to non- uniqueness for the Leray-Hopf solutions of the Cauchy problem for 
uq G L 2 . We plan to address these issues in future work. 

Our paper is organized as follows: in section 2, we prove an 'e- regularity' criteria for 
a generalized Navier-Stokes system; in section 3, we study the local in space near initial 
time smoothness of Leray solutions; in section 4 we study the asymptotics of forward self 
similar solutions to Stokes and Navier-Stokes equations; in section 5 we prove the existence 
of forward self similar solutions for large —1 homogeneous initial data. 

Notation: We use standard notations in our paper. For instance, Br{xq) denotes a ball 
centered at x with radius R in R 3 , Br : = Br(0); for z = (x , to), Q(R, ^o) : = Br(xq) X (t — 
R 2 , to], Qr := Q(R, (0, 0)); for any / in O, j a f := pr J f. We also use the following stan- 
dard notations in the literature: for vectors a and v, a <S> v is a matrix with (a ® v)y = a^-; 
for two matices a, b, (a : b) = a^b^ where we assume the usual Einstein summation con- 
vention; for a matrix valued function / = (fij), div / is a vector with (div /)$ = (^V djfij); 



(u)q(r, Z0 ) ■= f Q{Rzo) udz, [u) r := (u) Qr ; (p)b«(xo)(*) : = fB R(xo) P( x ^) dx ^ (p)r(*) := (p)bt(*); 
Y(u,p,Q(R,z Q )) := (f Q[RtZo) \u - (u) Q{RtZo) \ 3 dz) 1 ^ + R{j Q \ R>Zo) \p - (p)b k( , )(*)| 3/2 ^) 2/3 ; 



Y(u,p, Qr) := Y(u,p, Q(R, (0, 0))). We use C to denote an absolute and often large positive 
number, c a positive small absolute number, e the positive small numbers, C(a, /?,...) when 
the number depends on the parameters a, f3, .... Cp ar (C) denotes the Holder space with 
respect to the parabolic distance when O is a space time domain. We adopt the convention 
that nonessential constants can change from line to line. We use uq as a divergence free 
initial data throughout the paper, unless defined otherwise. 



Our goal in this section is to prove an e-regularity criteria similar to that of Caffarelli- 
Kohn-Nirenberg for a generalized Navier-Stokes equation. Our setting is as follows: 
Let O be an open subset of R% x R t) a G L^ C (C) with m > 5 (not necessarily an integer), 
div a = 0. We call a pair of functions (u,p) suitable weak solution to 



2. e- REGULARITY CRITERIA 



d t u — An + a ■ V« + div (a <S> u) + u ■ Vu + Vp = 

div u = 




(2.1) 
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if u G L™L 2 X n L 2 t H\{0') for any open subset O' C O' m O, p e L^(0), such that (u,p) 
satisfies equations ( 12.1 p in the sense of distributions in O, and 



\u\ 2 . \u\ 2 2 ■• /I""' 2 



m + a) J + m div (a (g) w) + div (wp) < 0, (2.2) 

in the sense of distributions. Recall that a distribution v in (9 is called nonnegative if 
(t> ; 0) > for any G C%°(0) with > 0; u div (a £g> m) is a distribution with 



(w div (a<8>w),0) — — diUjdjUi(J)(x,t)dxdt — / a,iUjUidj(j)(x,t)dxdt. 

Jo Jo 

The terms in 12.21 make sense due to the regularity assumptions and u G L^/. 3 (0) by known 
multiplicative inequalities. 

The main theorem in this section can be stated as the following: 
Theorem 2.1. (e- regularity criterion) 

Let (u,p) be a suitable weak solution to equations K2. 1\) in Q x with a G L m (Qi), m > 5, 
div a = 0. Then there exists eo = e (m) > with the following property: if 

/ r \ V 3 / /> \ 2/3 / » \ 1/m 

^y- |w| 3 etedtj + ^y |p| 3/2 dxcftj + yj \a\ m dxdtj < e , (2.3) 

then u is Holder continuous in Q1/2 exponent a = a(m) > and 

\\u\\ c ^ {Ql/2) <C(m,e ). (2.4) 

Remarks: The proof of this theorem follows the general line of presentation in [7tlT5|[TS]. 
There are some additional complications due to the new terms a ■ V« + div (a <8> u) as we 
shall see below. 

Before going into the proof of the theorem, we need the following two lemmas to be used 
below. 

Lemma 2.1. Let f be a nonnegative nondecreasing bounded function defined on [0, 1] with 
the following property: 

for any 3/4 < s < t < 1 and some positive constants < 6 < 1, M > 0, /3 > 0, we have 

f(s) < 0f(t) + j^y. (2.5) 

Then, 

sup f(s)<C(0,/3,M), (2.6) 
«e[o,3/4] 

for some postive constant depending only on 8, (3, M . 

Remarks: The lemma is well-known, one can find a proof for example in [8]. 
Our next lemma is an estimate of a generalized Stokes system. 
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Lemma 2.2. Let a G L' m (Qi), with div a = and (Jq 1 \a\ m dxdt) l l m < M , for some positive 
M > 0, m> 5, let A G iT ; |A| <M,f = (fa) G L m (Qi) m& (/ Qi |/| m dz<ft) 1/m < M. lei 
u G LfL 2 x n L 2 t Hl(Qi) and p G L 3/2 (Qi) wit/i 

^ |u| 3 dzd^ 7 + ^ |p| 3 / 2 dxdt^j 1 < M, (2.7) 



Assume (u,p) satisfies 

in Qt, (2. 



<9tn — Au + a • V« + A • V« + div (a ® w) + = div / 

div it =0 

m t/ie sense of distributions. Then u is Holder continuous in Q±/ 2 exponent a = a(m) > 
and 

|h|| c?ar(Ql/2) <C(m,M). (2.9) 

Proof: The proof is based on a standard application of boostraping arguments. We sketch 
some of the details below. In order to use the boostraping argument, we need to show, 
suppose u G L q (Qu), q > 3, then u G L^(Qr-s), for some q > q. Here we can assume 
R > 3/4 and 5 is a small positive number. Let us rewrite the equations of (u,p) as 

d t u — Au + Vp = div (/ — a®u — u® a — u® X) \ . „ 
divu =0 J mQn - 

mq 

By Holder inequality, we see h := / — a®u — u® a — u®\ G L m +i (Qr) ■ Taking divergence 
in the first equation, we obtain 

Ap = div div (/ — a ® u — u® a — u ® X) in Q R . 

Set 

Pi = A _1 div div ((/ — a®u — u® a — u® X)xb r ) , 

and write p = p\+p 2 . Then Ap 2 = in Qr. Recall that R is in [3/4, 1]. By elliptic estimates, 
we get 

llpill -2^2_ < C\\h\\ JHQ 
Since > 3/2, we see p 2 verifies estimate 

with S being a small positive number. Then, 

d t u - Au = -Vpi - Vp 2 + div h in Qrs/2, 

where pi, p 2 and h satisfy above estimates. For a smooth cutoff function r\ with rj = 1 in 
Qr-36/4: and n = outside Qrs/2, set 

«iM = f e A ^[-V(p lV )(; s) + div (h V )(-,s)}ds, 
J —00 

VzM = [ e A ^V(p 2V )(-,s)ds. 



00 
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Write u = V\ + v 2 + v 3 . By estimates of heat equation, we see ||u2||i«>(Q il _ s ) < C(S, M). As 
for vi, by Young's inequality and the properties of heat kernel, we see v± G L r (Qn-s) for any 
r > such that 

1111 
_>_+__ 

r q m 5 

Since v% satisfies heat equation in Qr-zs/^i we see v% is smooth in Qr-s- Thus in summary, 
we get, 

u G L^Q R ) implies u G L r (Q R ^ s ) with ± = ± - ±(| - £). 

Since m > 5, after applying this boostraping argument for finitely many times, we can 
conclude u G L r °(Q 5 / 8 ) with r sufficiently large such that 

\a\\u\EL^(Q 5/8 ). (2.10) 

Then we can go back to the decompositions vi, V2, and V3, it is not difficult to verify that 
all of them are Holder continuous in Q1/2 with exponent a = a(m). If we keep track of the 
constants in the above process, it's clear we also have the estimates claimed in the lemma. 
Alternatively, one can use the closed graph theorem with appropriate function spaces to 
obtain the estimates, we omit the details here. The lemma is proved. 

Now we can return to the proof of Theorem 12.11 We first prove the following 'oscillation 
lemma', which roughly speaking asserts that if u is of 'small oscillation' in Qi, then the 
oscillation is even smaller in Qg for 6 < 1. 



Lemma 2.3. (Oscillation lemma) 

Let (u,p) be a suitable weak solution to equations \2. 1\) in Qi with a G L m (Qi), m > 
5, div a = 0, ||a||z, m (Qi) < c, |(tt)i| < M, for some small absolute number c > 0, and 
some positive number M. Then for any 9 G (0,1/3), there exists an e = e(9,M,m) > 0, 
Ci(M, m) > 0, and a = a(m) > such that if 

y(u,p,Qi) + |(u)i| \J- \a\ m dxdtj <e, 

then 

Y(u,p,Q e ) < C 1 (M,m)9 a nrfapM + |(u)i| (J- \a\ m dxd?j ' 

Proof: Suppose the lemma is false. Then there exists (ui,pi) and ai with the following 
properties: 

\{ui)i\ < M, \\ai\\ L m^ Ql) < c, div a { = 0, 

Y(ui,Pi,Q\) + |(ut)i| I 7 \ai\ m dxdt\ = 6j — > + as i — > +00, 

Y(u i ,p i ,Q e )>C 1 (M,m)9 a e i , 
and (ui,pi) satisfies equations (12. ip and inequality (12. 2 p with a replaced by a^. 
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Set 



fi 



Pi ~ (Pi)i(*) 



Then we have = 0, div div fa = 0, 

/ /> \ 1/3 / /■ \ 2/3 /■ /. \ 1/m 

l-f lotted*) +l-f |gi| 3/2 cZa;dtJ + I -f |/i| m c£rcftj <l,and, 

v 1/3 / - n2/3 

\Vi - ( Vi )e\ 3 dxdtj + {J- \q t - (qMt)\ 3/2 dxdt ) > d(M, m)6 c 



Moreover, (vi,qi) satisfies: 

d t Vi - Avi + eiVi ■ Vvi + a { ■ Vvi + div (a* <g> u,) + div /; + (ui)i ■ Vvi + Vqi = V, , 

div Vi = y ' ' 

in the sense of distributions in Q\ and 

_L AUr (\^t(^.n,.^(n,\.^ n \\ 

(2.12) 



- + |V^| 2 + div (Hi (e^ + + a,)) 

+fj div (/j + an ® u») + div i;^ < 0, 



in the sense of distributions in Q\. Here again the terms make sense due to our regularity 
assumptions and the interpretation of Uj div (a^ £g> + /j) as the one below inequalities f 12 . 2 j) . 
Since u$ G LfL 2 D L 2 Hl(Qi) and i>j satisfies equations ( 12. lip , we can change the value of i>j 

on a set of measure zero such that £ — > is continuous from (—1,0) to L 2 (l?i(0)), the 

weak L 2 space. 

From inequality (12. 12j) we obtain, 

-(x, t)(j)(x, t)dx + / / | Vvi\ 2 (f){x, s)dxds 



(0) 2 J-UB 1 (0) 

1 <• w 2 



< / — — (d t (f) + A(f>)dxds + / / -— [(iti)i + a» + 6iVi\V(/)dxds 
J-iJb l (o) z J-iJb^o) 2 

+ / / [(/i + °« ® w «) : ( S ^ v i ( f ) + u i ® V0)]cfa;<is + / / qiViVcjydxds, 

J-iJbt(O) J-IJBt(O) 



for any > with e C c °° (^(0) x (-1,*]). 
Let us define 



Ei(r) = ess sup /" — ^— (x, t)cZx + / / | V^| 2 (x, s)dxds, (2.13) 
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for < r < 1. By known multiplicative inequalities we have 

IMIiio/» ( g p) < CEi{r). (2.14) 

Then for any 1/2 < T\ < r 2 < 1, if we choose nonnegative test function <fi with support in 
Q r2 appropriately, we obtain the following estimates, with the help of the above local energy 
estimates, Holder estimates, and the estimates on Vi,qi,af. 

E(n) < C r)2 + -^-f l ^f(M +\a i \+e i \v i \)dxds 

+ f \fi\\Vvi\ + \ai\\vi\\Vvi\dxdt -\ — — f \fi\\vi\ + \ai\\vi\ 2 + \qi\\vi\dxdt 



< 



x 1/2 

C{M) 



(r 2 -n) 2 



+ / \fi\ 2 dxdt\ E(r 2 ) 1/2 + \\a i \\ LHQr2) \\v i \\ L io/3 {Qr J\Vv\\ L 2 {Qr ^ 



< ( C{M \ 2 + (C\\a t \\ LmiQl) + l/2)E(r 2 ). 

Note that we have ||oi||L m (Qi) < c with c small. So if we choose c such that Cc < 1/2, then 

we can apply Lemma [2.11 and conclude that E(3/4) < C(M,m). That is, Vi are uniformly 

bounded in L^°L 2 n L 2 H^(Q 3 /4). Thus by known embedding theorems and the fact that 

satisfies equations (12. lip (which provide crucial regularity in t), we can choose a subsequence 

of Vi (which we still denote as v^, such that for some A G -R 3 , v G L 3 (Q 3 / A ), q G L 3 l 2 {Q 3 /4) 

and a, / G L m (Q 3 / i ) with div a = 0, we have 

Vi(-,t) -± v(-,t) weakly in L 2 (B 3/4 ) for every t G (-(|) 2 ,0), 

Vi ->■ u strongly in L 3 (Q 3/4 ), 

% g weakly in L 3 / 2 (Q 3/4 ), 

(Mj)i — > A, Qj — ^ a weakly in L m (Q 3 / 4: ), 

fi^fm L m (Q 3/4 ). 

Moreover, we have 

f \v\ 3 dxdt] +( f \q\ 3/2 dxdt] +( f (|/| + \a\) m dxdt) < C, 
and |A| < M. 

From equations (12. lip for (vi,qi), we see 

d t v - Av + A ■ Vv + div (a ® v + v ® a + /) + Vq = 



Zl} '"ft/4. (2.15) 



By Lemma 12.21 on generalized Stokes system, we see for some a = a(m) > 0, v is Holder 
continuous in Q\/ 2 with exponent a, with respect to parabolic distance. More precisely, 

\v(x x M) ~ v(x 2 M)\ < C(M,m) (|xi -x 2 \ + |ti - t 2 \ 1/2 ) a . 
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Since Vi —> v strongly in L 3 (Q 3 / 4 ), we see 

\ 1/3 

\vi - {v^e^dxdt J < C(M, m)6 a , 

for i sufficiently large. 

Note that from equations (12. lip we have 

— Aqi = div div (ejt>j ® v $ + v i <8> Oi + <8> u») • 

Let <fc = + g 2 , where 

q] = (-A) _1 div div ((e^ <g> v { + ® a, + a, ® fi)XB 3/4 ) , 

with X-B3/4 being the characteristic function on -B3/4. Since V{ strongly converges to v in 
L 3 (Q 4 / 3 ), we have qj — q~i strongly converges to in L 3//2 (Q 4 / 3 ), where 

<?i = (-A) _1 div div ((u <g> a; + a 4 <g> u)xb 3/4 J • 
Since w is bounded, we obtain by estimates of Riesz operators q\ G L m (Qi/2)- Thus 

\ 3/2 

|*|* ixit ) 

l/m 



/ r \ V" 1 

<0(j IftptfedtJ <C(M,m)^ 5/r 



Therefore, for z sufficiently large, we have 

\ 3/2 



0^ |g^| 3/2 dxrf^ <C(M J m)6 1 



5/m 

\ 2/3 



By definition, Ag 2 = in Q3/4 and ( ig 3/4 l^f | 3 ^ 2 dxdt) < C . Thus by elliptic estimates, 
we obtain, 



0(J q \q*-(q*Mt)\ 3/2 dxdt^ < C6^ j jVqK-M%s 5/12) d^ 

< ce( e- 1 ' 2 f [ \ q f 2 dxdt\ 



< ce 2 '\ 

Therefore, summarizing the above, we see 

(J" \Qi ~ (q l )e(t)\ 3/2 dxd?j ' < C(M,m)e min ^ 3 ' 1 ^ m K i 

for i sufficiently large. This, together with the estimates on Vi, shows 

Y( Vi , qi ,Qe)<C(M,m)e a , 
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for i sufficiently large, if we choose a{m) sufficiently small. This contradicts Y(vi,qi,Qg) > 
Ci(M,m)6 a if we choose Ci(M,m) > 2C(M,m). Thus the lemma is proved. 

The above lemma admits the following iterations. 
Lemma 2.4. (Iteration of the oscillation lemma) 

Let (u,p), M, e(9,M,m), Ci(M,m), a(m), c and a be as in the above lemma, with |(m)qJ < 
M/2. Let (3 = a/2. Choose 9 e (0,1/3) such that C 1 (M,m)9 a - fl < I, and 9 < a with 
c\ = ci(M,m) being some small number. Then there exists e*(9, M,m) sufficiently small, 
such that if 

/ r \ 1 / m 

Y(u,p, Qi) + M yj- \a\ m dxdtj < e*, (2.16) 

then for any k = 1, 2, . . . , we have 

KtOg^l < M, (2.17) 

(\ 1/m 
-f \a\ m dxdt) 9 k ~ l < e* < e{9,M,m), (2.18) 

Y(u,p,Q ek ) < 9 P ^Y(u,p,Q ek -i) + \(u) gk -i\ (j- \a\ m dxd?j 9 k ~ 1 \ . (2.19) 
Proof: We prove the lemma by induction. 

For k — 1, the conclusion follows from Lemma 12.31 if we choose e* such that e* < e(9, M,p). 
Suppose the conclusion is true for k < k , k Q > 1, we show it remains true for k = k + 1. 
By induction 

|(u)o, fc _xl <M, 

Y(u,p,Q 6 k-i) + \ (u)gk-i\ j- \a\ m dxdt 9 k 1 < e*, 

Y(u,p,Q ek ) < 9? (^Y{u,p,Q ek -,) + | | ^ [a] m tfa;dt^ 0*"^ < ^e*, 

for all k < k . Thus, 

Y(u,p,Qg k ) < 9P\Y(u,p,Qgk-i) + 9 k ~ 1 M(j' \a\ m dxdt 
\ \JQ gk -i 

( f r \ V"^ 

< 9^ \ Y(u,p,Q ek -i) + 9 {k ~ 1){1 - 5/m) M yj- \a\ m dxdt 

< 9^Y(u,p,Qg k -i) + 9 k ^e, 
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for all k < ho, with fli = min{/3, 1 — 5/m}. Simple calculations with a repeated use of the 
above inequalities show: 

Y(u,p,Q ek ) < 9 k ^Y(u,p,Q 1 ) + k9 k/Sl e„ Wk < k . 

Thus, 

ko 
k=l 

fco / , \V3 



fed 



< r 5 / 3 ^F( M , P ,g, fc - 1 ) + |( M ) Ql | 

fc=i 

< r 5 / 3 ^ (fl^^e* + e,(Jfe - l)^^- 1 ^) + M/2 

fc=i 

< r 5 / 3 (I - fl^e* + r^C^, 0) + M/2. 
If we choose e* = e*(#, M, m) to be sufficiently small, we see 

l(«)Q efeo l<M. 

Moreover, 

l/m 



Q^o) + O ko \(u) e k \ f \a\ m dxdt 
\JQ 9 k 

< 9 p e* + 0(i-6/«O*o e , < e*, 
if we choose 6 1 < c(M, m) to be sufficiently small. Set 



p(x,t) = -^rrqi^T-, -Jtrr), and 

I \ 1 7/ X * * 

°0M) = T^te, 
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One can verify that (v, q) is a suitable weak solution to equations ( 12.1 \ with a replaced by b 
in Q\. Moreover, 

= ^° ^>,P, Q fl *o) + ko \(u) QekQ | ^£ |ardxd^ j < e„, 

f r \ l i m Bk f c \ 

(j \b\ m dxdtj <6 ko -°-^\^J \a\ m dxdtj < c. 



Thus, by Lemma 12.31 we obtain, 



Y(v,q,Q g )<9^lY(v,q,Q 1 ) + \(v) Ql \ (/ \b\ m dxdt) ), (2.20) 



Qi / 

that is, 

Y(u,p,Q 0ko+ i) < 9? ^Y(u,p,Q ek0 ) + \(u) ek0 \ i^j \a\ m dxd^j 9 k ° \ . (2.21] 

The lemma is then proved. 

By translation and dilation, we obtain the following corollary. 

Corollary 2.1. Let (u,p) be a suitable weak solution to equations pL 1\) in Q(R,zq), with 
a G L m (Q(R, Zq)) , div a = 0, \(u)q^ Zo )\R < M/2, 9, (3 are as in the above. Then there 



exists e* = e*(0, M, m) such that 

RY(u,p,Q(R,z )) + RM (i 

\Jq 



l/m 

\a\ m dxdt I < e* 



Q(R,z ) 



implies, for k > 1 : 



R\{u)q{6^r,z )\ < M, and 
Y(u,p,Q(9 k R,z )) 

<9P{Y (u, p, Q(9 k ~ 1 R, *„)) + R9 k ~ 1 \ (u) Q . lR \ (1 \a\ m dxdt] 

\ \JQ(e k ' 1 R,zo) J 

Proof of Theorem I2.lt 

It is clear if we choose eo sufficiently small, we can apply Corollary 12. II in Q(l/2, zq) for any 
z o £ Qi/2- Note that \(u)Q gkR \ is bounded and m > 5. Thus we can conclude 

Y(u,p,Q(z ,Qg k )) < C(9,M,m)9 ka , 
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for some a = a(m), where we can choose M < 1, 6 — 6(M,m) = 9{m). (There is a slight 
abuse of notation, in particular, this a is smaller than those appearing in the oscillation 
lemma). In particular, 

\ 1/3 

\u - (u) Q{ek<zo) \ 3 dxdt < C(0, M, m)9 ka , 
'Q(e k ,z ) J 

for all zq G Q\/2 an d k > 1. By Campanato's lemma, we conclude u is Holder continuous in 
Qi/2- The theorem is proved. 



(/ 



In applications, it is cumbersome to have the "smallness condition" on a. We can remove 
this condition and get the following theorem. 

Theorem 2.2. (Improved e-regularity criteria) 

Let (u,p) be a suitable weak solution to equations Ii2. 1\) in Qi, with a G L m (Qi), div a = 0, 
ll a IU m (Qi) — M, for some M > and m > 5. Then there exists ei = ei(m,M) > with the 
following properties: if 

\ 1/3 / r \ 2/3 



C \ i f r \ 1 

f \u\ 3 dxdt I + I -r \p\ 3/2 dxdt J < ey 
JQi J \JQi J 



then u is Holder continuous in Q1/2 with exponent a = a(m) > and 

IMIc^Wy) < C ( m > ei, M) = C(m, M). (2.22) 

Proof: Choose < Rq < 1/2, a small positive number to be determined below. For any 
z o — ( x Oito) G Q1/2, we would like to apply a scaled version of Theorem 12.11 for (u,p) in 
Q(R,z ). Set 

u[x,t) = —v\ 

t \ 1 
POM) = -B2?( 

a(x,t) = —b(- 



x 


- x 


t 


-t 




Ro 


J 


R o 


X 


- x 


t 


— 1 , 




Ro 




Rl 


X 




t 


— to. 



Ro Ro Rq 

We see that (v, q) is a suitable weak solution to equations (12. ip with a replaced by b in Qi- 
Moreover, 

and 



^ D l— 5/mii 11 s~i r>l—5/ m n /r 

< K |M|L™(<2(i?o,zo)) ^ CK M > 



= Ro(i \u\ 3 dxdt] +(i \p\ 3/2 dxdt] R 2 

\Jq(Ro,z ) J \Jq(R ,z ) J 



'Q(R ,z ) / \JQ(R ,z ) 

< C(RqR 5 ^ 3 + R^Rq 10 ^ 3 )ei < Ci? 4//3 ei. 
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Thus, 

f \ f r \ 2 ^ 3 / f \ 

f \v\ 3 dxdt +lf \q\ 3/2 dxdt +lf \b\ m dxdt 
JQi J \JQi J \JQi J 

<Rl~ 5/m M + CR~ 4/3 e 1 . 

Thus, if we choose R such that R l Q ~^^ m M < e /2, fix R , R = R (M,m) and choose e\ 
such that CR~ i/3 e l < f . Then we can apply Theorem 12. II to (v, q) and conclude v is Holder 
continuous in Qi/2- Scale back and collect all constants, the theorem is then proved. 



3. Local in space near initial time smoothness of Leray solutions 

In this section, we use the 'e- regularity' theorem proved in the last section to study the 
local in space near initial time smoothness of the so called Leray solutions. Our setting is as 
follows. 

Let u G L 2 oc (R 3 ) with div u = and sup X0gJ?3 J Bl ( X0 ^ \uo\ 2 dx < oo. We recall the definition 
of Leray solutions in [19], see also [T2] . 



Definition 3.1. (Leray solution) A vector field u G L 2 0C (R 3 x [0, oo)) is called a Leray 
solution to Navier-Stokes equations with initial data Uq if it satisfies: 

i) esssup < t<fl2 sup X0&R3 J Br(xo) ^-{x,t)dx + sup^^s ff J Br(xo) \Vu\ 2 dxdt < oo, and 

lim / / 

\x \^°°J Jb r (x ) 

for any R < oo. 

ii) for some distribution p in R 3 x (0, oo), (u,p) verifies Navier Stokes equations 

d t u — An + u ■ Vu + Vp = 
div u =0 



-R 2 r 

lim / / \u\ 2 (x,t)dxdt = 0, 



in R 3 x (0,oo), (3.L 



in the sense of distributions and for any compact set K C R 3 , lim 4 ^ 0+ — u \\ L 2^ = 0. 

iii) u is suitable in the sense of Caffarelli-Kohn-Nirenberg, more precisely, the following local 
energy inequality holds: 

r°° r r°° r | M | 2 \ u \z 

/ / \X7u\ 2 (j)(x,t)dxdt < / —z~{dt4> + A0) + —p—u ■ V0 + pu ■ W<pdxdt (3.2) 
Jo Jr 3 Jo Jr 3 2 2 

for any smooth > with supp d R 3 x (0, oo). The set of all Leray solutions starting 
from uq will be denoted as N{uq). 

Remarks: For general existence result of Leray solutions, see [HUH [19]. ^or us ' ^ ne a 
priori estimates of Leray solutions below are more important, since in our situation when 
Mo is usually better than that in [19], the existence can be proved in simpler ways. In 
the case the initial data is in L 2 (R 3 ), the notion of Leray- Hopf weak solutions is often 
used (see [16] for example). The difference is that Leray-Hopf weak solutions belong to 
LfL 2 , n L 2 H X (R 3 x [0, oo)). It is clear that our definition includes such solutions. Note 
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that we impose a decay condition on u in i). This condition allows us to calculate p in 
the following way: \/B r (xo) x (0, £*) C R 3 x (0, oo), take a smooth cutoff function with 
4>\b 2 t(xo) = 1) then there exists a function depending only on x ,r,t,<p (we suppress the 
dependence on xo,r, in our notation) such that for (x,t) G B r (xo) x (0, i*) 

i) = — A _1 div div(w ® «</>) — / {k(x — y) — k(x — y)) u <g> u{y, t) (1 — 4>(y)) dy + p(t) 

(3.3) 

where k(x) is the kernel of A x div div. 

The right hand side is well defined since u satisfies the estimates in i) and 

\k(x - y) - k(x -y)\ = 0(- — tt) as \y\ -> oo. (3.4) 

Fo - y\ 

The situation is similar to extending the domain of singular integrals to bounded functions, 
see for example [19] and [23] . 

For Leray solution u G Af(uo), we have the following a priori estimates, first proved in [T9] . 
see also a simpler proof in [12]. These estimates have played an important role in [1~2||22]. 
see also 1251. 



Lemma 3.1. (A priori estimate for Leray solutions) 

Let a = sup XogR 3 f BR r XQ \ ^-(x)dx < oo for some R > and let u be a Leray solution with 
initial data u Q . Then there exists some small absolute number c > such that for A satisfying 
< A < cmin{ct" 2 i? 2 , 1}, we have 

f \ u ? f XR2 f 

ess sup sup / —— (x, t)dx + sup / / | Vu| 2 (x, t)dxdt < Cot. (3.5) 

0<t<AR 2 xoS-R 3 Jb r (x ) 2 x S-R 3 JO J B R (x ) 

Remarks: Note that from the formula (13. 3p and the a priori estimate of u, we get the 
following estimate for p which will be useful: 

sup / / \p-p(t)\ 3/2 dxdt<Ca 3/2 R l/2 . (3.6) 

x &R 3 JO J B R (x ) 

In the above estimate on p, more precisely, p(t) = p Xo ,R(t). That is, we need to choose some 
appropriate constants p Xo ,R(t) to satisfy the inequality. The point here is that such constants 
depending on xq, R, t exist. This remark is effective throughout the paper. 

Now we can prove our our first important result. 



Theorem 3.1. Let u G Lf oc (R 3 ) with sup Xo£R s f Bl r XQ \ \u\ 2 (x)dx < a < oo. Suppose u is in 
L m (B 2 (0)) with ||^o||L m (B 2 (o)) < M < oo and m > 3. Let us decompos^ u = + u 2 , with 
div = 0, u\\b a/3 = uo, supp u\ d B 2 (0) and ||itolU m (ii 3 ) < C(M,m). Let a be the locally 

1 Such decomposition is well-known. One can for example first localize uq using a smooth cutoff function, 
and then use Bogovskii's lemma to deal with the divergence-free condition. See for example [IJ[9]. 
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in time defined mild solution to Navier-Stokes equations with initial data u\. Then there 
exists a positive T = T(a, m, M) > 0, such that any Leray solution u G N{uq) satisfies: 
u- a G C^ T (B 1/2 x [0,7]), and \\u - a|| c ? ar (5^x[o t]) < C(M,m,a), for some 7 = 7(771) G 
(0,1). 

Remark: We can certainly choose T(M) > such that a is defined on R 3 x [0,T(M)]. The 
point of the theorem is that regularity of solution to Navier-Stokes equations depends locally 
on initial data, as least when Holder continuity is concerned. 

Proof: By assumption a solves the Cauchy problem for Navier-Stokes equations with initial 
data Uq in R 3 x [0,Ti], where T\ = Ti(M,m), namely: 

d t a - Aa + a ■ Va + Vp =0 



dJf^o} i-^xM), (3.7) 



and a(-,0) = itj . (3.8) 

It is well-known how to construct the so called mild solution to Navier-Stokes equations, see 
for example [T3"l[To'll20"] . In our case, it is even simpler, since Mq G L m with m > 3 is subcrical 
with respect to the natural scaling of the equation. We can follow the arguments in the 



Appendix of 171, and obtain a G V"^ (R 3 x (0, Ti))) with llall 5m, , , < CM. Note that 

> 5 since m > 3. Moreover, by the estimates on a and by treating the nonlinear term as 
pertubation, we can recover a local energy estimate for a: 

sup f — — (x, t)dx + f f \Va\ 2 (x,t)dxdt < C(M,m), 

<t<T 1 JB 1 (x ) 2 Jb^x^JO 

for any x G R 3 . Write u = a + v , we can verify that v satisfies: 

d t v — Av + v ■ Vf + a ■ Vv + div (a ® v) + Vg = 1 , , 

divv = J ^' 9j 

in the sense of distributions in R 3 x (0, Ti), here q = p — p with p being the associated 
pressure for u; and the local energy inequality 



ess 

0<t<Tj 



M 2 









71 

d t ^ A^- + I Vw ' 2 + div 2 ^ + a )) + u div ( a ® w ) + div < °> 

in the sense of distributions in R 3 x (0,Ti); 

Note also that ^ol-B 4/3 = 0. Since (it,p) satisfies the a priori estimates in Lemma I3TT1 (and the 
remarks below it), (a,p) is regular, we obtain the following estimates for (v,q) in ^(0) x 
[0,T 2 ), T 2 = T 2 (a,M,m): 

ess sup - / \v | 2 (x, t)dx + / / | Vw| 2 (x, s)dxds 
o<t<T 2 2 Jb 2 (o) Jo Jb 2 (o) 

+ (/ / \q\ 3/2 dxds} <C(a,m,M). 
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From the local energy inequality for v, and lim^ 0+ t) ||i 2 (B 4 y 3 (o)) = 0; we obtain 

1 <• 



v\ (x, t)(J)(x)dx + / / |Vf| (x,t)(j)(x)dxds 



b 4/ / 3 J J -B4/3 



2 

< I I ^A(j)dxds+ I I ^-(v + a)V(j>dxds 

Jo J B 4/3 2 Jo J B 4/3 2 

+ f I [a <g> u : (Vf0 + u <g> V0)] + <?f • V(j)dxds, 
Jo Jb a/3 

where G C?(B i/3 ), <f>\ Bl = 1, > 0. 
By multiplicative inequalities, we know 



T 2 f x3/10 

| 10/3 (M <C(a,m,M). 



'B 2 (0) / 

Thus from the above, we see by Schwartz inequality: 

\( \v\ 2 (x,t)dx + [ [ \Vv\ 2 {x,s)dxds<C{a,m,M)t rain{1/30 ' n ^ } , 
2 J_Bi(o) Jo ifii(o) 

for t < T 2 . From 

Ag = — div div (t> ® v + a (g> t> + v <8 a), 

we can see q G -£^/ c 3 . Thus 



loc 



< C(a,m, M)t 1/15 . 

The importance of these estimates lies in the fact that they provide crucial "quantitative" 
information on the decay in time as t — » 0+. Now for to fixed, whose precise value is 
to be determined later, extend v, q to -Bi(O) x (—1 + to, to] by setting v = 0, q = for 
(x,t) G B 1 x (-1 + t ,0]. Extend a to #i(0) x (-1 +t ,^o] by setting a(t,x) = for t < 0. 
The extended function (w, g) is a suitable weak solution to the generalized Navier-Stokes 
equations (12. ip with the extended a in -Bi(O) x [— 1 + t ,to\- Note here that 

hm_ \\v(-,t)\\ L 2 [Bim = 

plays a crucial role: it guarantees that dtv and df^- will not cause any problem across t — 0. 
Then clearly if we choose £ = ^o( a ; m i M) sufficiently small, we can apply Theorem 12.21 and 
conclude v is Holder continuous in B\/ 2 x [0,t ], with ||w||c5'ar(Bi/2x[o ) to]) — C(a,m, M), for 
some 7 = 7(771). The theorem is proved. 

For applications below, we state the following simple (and certainly well-known) lemma 
for heat equation without proof. 

Lemma 3.2. We have the following estimates: 

1. If u G CP(R 3 ) for some f3 G (0, 1), then e At u (x) G C^ r (R 3 x [0, 1]), with 

\\e At u (x)\\ c ^ {nixm) < C\\u \\ c p m . (3.10) 
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2. If f E L°°(R 3 x [0, 1]), then f* Ve A(t_s) /(-, s)ds E C^ T (R 3 x [0, 1]) for any (3 E (0, 1), and 
|| / Ve^f(-,s)ds\\ 

Cg^R? x[0,l]) ~ C (P)\\f\\L°°(R3x[0,l])- (3-11) 

J 

The above theorem implies the following result. 
Theorem 3.2. (Local Holder regularity of Leray solutions) 

Let u E Lf oc (R 3 ) with sup x . ogR3 f Bl r X0 ) \u\ 2 (x)dx < a < oo. Suppose u is in C 7 (I? 2 (0)) with 
|l M o||c^(s 2 (o)) < M < oo. Then there exists a positive T = T(a, 7 , M) > 0, such that any 
Leray solution u E M{uq) satisfies: 

u E q m{ (FhJl x [0,T]), and |k|| CpV(5 ^ x[0iT]) < C(M,«, 7 ). (3.12) 

Proof: Let us decompose uo = u\ + Mq with div u\ = 0, Mo|b 4/3 (o) = uo, supp u\ <s ^(0) 
and 1 1 Uq \\c-y(R?) < CM. Let a be the mild solution to Navier-Stokes equations with initial 
data Uq in R 3 x (0,T(M)). Then Theorem 13. II implies that u — a is Holder continuous with 
some exponent j3 E (0,7) in B 1 / 2 x [0,T] with some T = T(a, 7, M) G (0,T(M)). Since the 
initial data Mq for a is in C 7 , it is not difficult to show that a E C 7 (R 3 x (0,T)). Thus u 
is Holder continuous with exponent j3 in Bi/ 2 x [0,T(M)] . By using a routine boostraping 
argument, one can improve the exponent to 7. Since this argument will be used one more 
time below, we sketch some of the details here for the reader's convenience. Note that u 
is Holder continuous in Bin x [0, T], thus from the representation formula (13. 3p for p and 
estimates for Riesz transform, we know p is bounded in -B7/16 x [0, T] modulo some function 
p{t). Now rewrite the equation for u as 

d t u — Au = —div (u (2) u) — Vp. (3.13) 

Choose a smooth cutoff function rj with 77 = 1 on B 3 /g and rj = outside -67/16- Write 

ui{-,t) = [ e A( '^ s) [-div (u®wn) - V{pr])}{-, s)ds, 
Jo 

U2(-,t) = e At (u r]). 

Let u = ui + u 2 + M3. By Lemma I'd. 21 we see u\ and u 2 are Holder continuous with exponent 
(3. Note that u 3 satisfies 

d t u 3 -Au 3 = in B 3/8 x [0,T], 

and «3(-, 0)|b 3/8 = 0. Thus M3 is smooth in 5i/ 4 x [0,T]. In summary tt is Holder continuous 
in B1/4 x [0, T] with exponent /3. Then the theorem is proved. 

4. Estimates of forward self similar solutions to Navier-Stokes and 

Stokes equations 

In this section, we start to study forward self similar solutions to Navier-Stokes equations 
and a related nonhomogeneous Stokes system. Our setting is as follows. 
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Let u be a Leray solution with initial data u . Suppose Xu (Xx) = u (x), Xu(Xx,X 2 t) 
u(x,t) for any A > 0. We also assume Wo|asi(o) £ C°°(dBi(0)). Then it is easy to see 

|V"«.h)l < V |a| > 0. 



\x 



Our first main result in this section is the following theorem. 
Theorem 4.1. (A-priori estimate for forward self similar solutions) 

Let u, uq be as in the above. Then U(-) := u(-, 1), the solution profile at time t = 1, belongs 
to C°°(R 3 ) and 



C(a,u ) 



|0° (U(x)-e*u (x))\< {i ^^ H , V 



H>o. (4.r 



Remarks: Here and below, constants C(u , . . . ),T(uq, ...)... only depend on the magni- 
tude of u Q and its finitely many derivatives on the unit sphere. Similar estimates with more 
precise asymptotics have been proved in [2] when the initial data is small in appropriate sense. 



Proof: Apply Lemma [3. II with R = 1, we see (set M := ||wo||c(SBi 



sup -f \u(x,t)\ 2 dx+ [ [ \Vu(x,t)\ 2 dxdt<C(M), T 1 = T 1 {M). (4.2) 
o<t<Tj 2 Jbi(o) Jo Jb!{o) 

For fixed t* < T±, with to be determined later, since u(x,t) = ^w(^=, 1) = 77J ^ ( "3* ) ' we 
have 

C(M) > 1/2 f \u(x,Q\ 2 dx + [ [ \Vu{x,t)\ 2 dxdt 



> [ \u(x,l)\ 2 dx + y^ I 



B i (0) ° JB^_(0) 

•JtZ 



> f \U(x)\ 2 dx + ^ [ \VU(x)\ 2 dx. (4.4) 

2 Jb^_(0) ° ■/ -B i (0) 

On the other hand, for Vxo, |a?o| = 8, since u G C^-B^o)), we can apply Theorem 13. II and 
some simple boostraping arguments to show the following: 
there exists T 2 = T 2 (M) > such that V a, 

\\ d tdy\\ L °°(B 1/8 {x )xio,T 2 ]) <C(a,uo), 
this is true for any u G M{uq) . 

Since V A > 0, Aw(Ax, X 2 t) is also a Leray solution with initial data Uq, we obtain 

\\\ a \+ 1 d a u{Xx , XH) - d a u {x )\ < C{a, u )t, 
for any A > 0, |a| > 0, t < T 2 (u ). 

Take A = we obtain \{±)\ a \+ l d a u{^ 1) - d a u(x ) \ < C(a,u )t. 
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Setting y — and using the homogeneity of d a u , we get 

\d"U(y)-d"u (y)\<^^, V\y\>^=. (4.5) 

\y\ v +2 

Now choose sufficiently small, t* = £*(M), we see from inequality (14. 3ft : 

(\U(y)\ 2 + \VU(y)\ 2 )dy<C(M). 

Since u(x,t) satisfies Navier-Stokes equations, it is easy to verify U satisfies 

-AU-2-VU-% + U-VU + VP =0 

div U =0 

Thus elliptic estimates give 



B 

2 



in R 3 . (4.6) 



\M-)\\cHB , )<C{k,M)- 
These estimates, combined with the properties of heat equation, finish the proof. 

For later use, let us study a nonhomogeneous Stokes system with singular forcing. Our 
result is the following Lemma. 

Lemma 4.1. (Decay for the linear singularly forced Stokes system) 

Let f £ C(R 3 ), suppose v £ LfL1(R 3 x (0, T)) for any T < oo, and some 7 > 1, suppose v 
satisfies 

*,-a„ + v„ =r^)i mR3x{0:OO)> (47) 

div u = J 

/or some distribution p, and liim_ s>0 + \\v(-, 0IUt(r 3 ) = 0. T/ien 

i) i/0 also satisfies the above conditions, then v — v. 

ii) if f satisfies M := sup ;EgR 3(l + |x|) 3 |/(x)| < 00, then 

«(-,*)= j\ A ^P-Lf(-=)ds, (4.8) 

where P is the Helmholtz projection operator. Let V(x) = v(x, 1), then \\V\\ci,<*(b r ) — 
C(a,R)M for a £ (0,1) and 

sup ((1 + \x\) 2 \V(x)\ + (1 + |a;|) 3 |W(x)|) < CM. (4.9) 

Hi) if f satisfies M : = sup x6i? 3(l + |x|) 4 |/(a;)| < 00, then 

v(-,t)= f e^P^rJi^ds. (4.10) 
Jo s 1 V s 

Let V(x) = v(x, 1), then HV^Ici^^) < C(a, R)M for a £ (0, 1) and 

sup ((1 + \x\f\V(x)\ + (1 + \x\Y\VV(x)\) < CM. (4.11) 

xeR 3 
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Proof: The uniqueness is easy. We only need to show that if / = and for some 7i, 72 > 1, 

v G (LI 1 + Ll 2 )(R 3 x (0,T)) for any T > and, 
t Um \\v(t, OH (i#+i#)m = 0, 

then v = 0. Set u = curl u, then d t u — Aoj = in R 3 x (0, 00). Since 



ttm ||-y(-,t)||(L71 + L72)(fl3) = 0, 

we can extend u to R 3 x i? by setting u; = for v < 0, and the extended function, which 
we still denote as u, satisfies dtoo — Au = in R 3 x R. Here again there is no problem 
showing that the equation is satisfied across t = since u decays to as t — > 0+. One can 
for example first mollify u in x and the mollified function is smooth in both x and t. For 
the mollified function the claim is clear, then we can pass to the limit to show our claim. 
Since we have bounds for u in some negative Sobolev space and uj = for t < 0, we conclude 
u = 0. Thus Av = in R 3 x (0, 00). Therefore v = 0. 

Let us now prove part ii) and part iii). By the uniqueness result, we only need to prove the 
claimed estimates. Denote the kernel of Pe A by k(x), then k(-) G L 1+e (R 3 ) for any e > 0. 
By Young's inequality it is easy to get 

"t 

l+e 



< C(e) f \\{t - S )^ 3/2 A;(-=)|| Li+ ,|| S - 3 / 2 /(^=)ll^^ 



< C(e)Mt l - 



36 

l+e 



Thus, 



Now let us prove the decay estimates of V. The proof is a direct consequence of the following 
inequality (which can be proved by simple calculations) with a, (3 = 3, 4 and R := \x\ > 8 : 

lf 1 1 dvdt<l R ~ 3l ° gR if « = / ? = 3 ' (4121 



'0 J& (\x -y\ + vT=t)<* (\y\ +Viy "I R~ a ~ p+i otherwise 
For part i), we have 

1^)1 * HL j^hr^w mhw ivit - w_3logW 
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for \x\ > 8. For part ii), we have 

for \x\ > 8. Thus the decay estimates are proved. Since V also satisfies an elliptic equation: 



-AV-f -W-^ + VP =f 



2 div V = 



in P 3 , (4.13) 



the estimates in Bc(0) is simple. 



5. Existence of forward self similar solution for large initial data 

In this section, we prove the following theorem. 

Theorem 5.1. Let u G C°°(P 3 \ {0}) satisfy Xu (Xx) = u (x) for all A > 0, div u = 0. 
Then there exists u G C°°(P 3 x (0, oo)), with Xu(Xx,X 2 t) = u(x,t) for all A > 0, and 
u G M{uq), that is, u satisfies 

d t u - Au + u ■ Vu + Vp = 1 . „ 3 , w , . 

jj v u — o f m ' ' °°' * or some P- (5-1) 

Moreover, let U(x) = u(x, 1), then 

Iff- (i/W - «%,(*)) I < (^'ffw 

Proof: By Theorem 14.11 it suffices to show there exists u G M{uq) with the scaling 

Au(Ax, A 2 t) = u(x, t) for all A > 0. (5.2) 

Denote 

X = {V G (^(P 3 ) : div V = 0, sup ((1 + \x\) 2 \V(x)\ + (1 + |x|) 3 | W(x)|) < oo}. (5.3) 

For any V G X, we define a natural norm 

\\V\\ X = sup ((1 + \x\) 2 \V(x)\ + (1 + |x|) 3 |W(x)|) . (5.4) 

Set C/q = e A UQ. Introduce a parameter /i G [0,1], set £/o M = hUq. We will follow Leray's 
method to prove the existence of u G J\f(u ) with Xu(Xx, \ 2 t) = u(x, t) for all A > 0. Due to 
the scaling invariance of u(x, t), we are essentially seeking the profile function U (x) = u(x, 1), 
where U(x) satisfies 

u - f • W 

div U =0 



-AU + U - VU - I - f • W + VP =0 



in P 3 , (5.5) 



and the correct asymptotics at spatial infinity. We will solve U in the following form 

U = U 0ll + V, where V G X (5.6) 
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It is clear u(x,t) = ^P(^) G J\f(fiu ) if and only if U(x) satisfies the above elliptic system 
and U (x) = Uq^ + V for some V G X, by Theorem 14.11 Thus we have reduced the problem 
to to finding V G X, with 

-Ay + v ■ w + u 0li ■ w + v ■ vp 0At - \ - § ■ w + vp = -p 0m • vp 0m 

div y =0 

in P 3 . We rewrite the above as: 

- AV - ^ - | ■ W + VP = -V ■ W - C/ „ ■ W - V ■ VU , - U , ■ VUo M - (5-8) 
Since 7 V satisfies the above equation if and only if v(x,t) := ^f^(^) satisfies 

divv =0 > (5.9) 

<,0) =0 J 

where 

P — ~V ' — ■ VV — 1/ ■ VPqm - Pom ■ VP 0m (5.10) 
has the decay properties in Lemma |4~T1 Thus for such P, equation (15. 9p is uniquely solvable, 
we denote the solution profile at time 1 as G(F) G X. This enables us to consider the 
following equivalent formulation, 

find V G X with V = Q(-V -VV -U^-VV -V ■ VP 0m - P 0m • VP 0m ). (5.11) 

Let K : X x [0, 1] — » X be defined as: 

V V G X, /x G [0, 1], tf(V, /i) := 6?(P 0M VP 0At ) + 6?(Po M VV + F VP 0M + VVV). (5.12) 

The first term on the right hand side is one dimensional. The second term by estimates in 
Lemma \4. II is compact. The compactness is due to the local C 1,a estimates and the fast de- 
cay at inifinity. Thus we conclude K G C 1 (X x [0, 1]) is compact. Therefore we are reduced 
to solve the following abstract problem: 

find V G X, such that V + K (V, fi) = 0, where \i G [0, 1]. 

At this stage, we are in a position to apply Leray's method, see for example [2T]. We need 
the following conditions to be verified: 

1. Solvability for \i small. This is already done, for example in j5j[TT] , note that it also 
follows from a simple implicit function theorem in our formulation. In the language of Leray 
Schauder degree theory, we can verify d(I + K(-, /x), 5^(0), 0) = 1 for \i small and some 
fixed M > 0. 

2. A priori estimate for solutions. This is done, in Theorem 14.11 

3. Compactness and continuity of K. This follows from the estimates of Q. 

Thus we can apply Leray's method, and conclude that for each \i G [0,1], there exists a 
solution V G X to V + K(V, //) = 0. Take /i — 1, the theorem is proved. 

With the existence theorem for smooth (away from 0) —1 homogeneous initial data, we 
can obtain existence results for not so smooth initial data. We illustrate the method with 
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Holder continuous (away from 0) initial data, although more general initial data can be 
considered. 

Theorem 5.2. Let u £ C" oc (R 3 \ {0}) with a £ (0, 1), Xu (Xx) = u (x) for all X > 0, and 

div u = in R 3 . Denote M = ||Mo||c* a (asi)- Then there exists u £ Af(u ), and u satisfies 
u{x,t) = Aw(Ax, XH) for all A > 0. Moreover, let U(x) = u(x, 1). Then U £ C°°(R 3 ) with 

\U(x)-e*u (x)\< {i ^ i+a . (5.13) 

Proof: Let us choose u% £ C°°(R 3 \ {0}) with Xu e (Xx) = u e (x) for all A > 0, div u e = 
in R 3 , ||wollc a (aBi(o)) < CM, and ||mq — WoU^aex) — > as e — > 0+. We can construct 
such Uq by first mollifying no on the unit sphere and then using the scaling invariance and 
applying Helmholtz projection operator to form Uq. We only note that the scaling invariance 
is preserved by the Helmholtz projection. By Theorem I5.1[ we can find u e £ N"(uq) with 
Xu e (Xx, XH) = u e (x, t), for all A > 0. Let U e (x) = u e (x, 1), then u e (x, t) = -L[/ £ (^). For any 

xo £ R 3 with |xo| = 8, since w £ C^^B^xq)) with ||moIIc q (b 4 (xo)) < C(M), by Theorem 13. 2\ 
there exists T(M) > 0, such that u e £ C^ r (B 1/2 x [0,T(M)]) and ||w £ ||c£ ai .(B 1/2 x[o,T(A/)]) < 
C(M). Thus, 

|^f/ e (^) - u e (x )\ < C{M)t a ' 2 , for t < T(M). (5.14) 
yt yt 

By the homogeneity of u e , we get 

\U e &) - u e Q &)\ < C{M)t 1 / 2+a /\ for t < T{M). (5.15) 
yt yt 

Notice that |xo| = 8 is arbitrary, we get 

\U e {x) - u%{x)\ < for |x| > Ci(M). (5.16) 



Moreover, by following the same arguments in the proof of Theorem 14.11 we can obtain 

\\U% ck{BRm < C(k,M,R) for Vi? > 0. (5.17) 
By combining the above estimates and using elementary properties of heat equation, we get 

\U e (x) - e A u e (x)\ < , °[ M l^ , for x £ R 3 . (5.18) 

Note also that since u e satisfies the Navier-Stokes equations for t > 0, U e satisfies 

-AU e + U e ■ VU e - f • W £ - + VP e = 1 . d3 

divtf' =0 j mR - 

By the estimates on U e , we can pass to a subsequence — > 0+, such that U £i — > U in 
C 2 (5 fl (0)) for all .R > 0. Thus C/ satisfies 



-At/ + £7 • VZ7 - f • V£7 - f + VP =0 



div C/ =0 



in R 3 , (5.19) 
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and 

\U(x) - e A u (x)\ < , ^f^L i for all x G R 3 . (5.20) 
ivy uv ;i - ^ + i^^i+lal v ; 

Setting u(x,t) = 77^(7%) j we can easily verify that u satisfies all the conditions in our the- 



orem. 
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